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Ca^^-activated Cl~ channels (CaCCs) are involved in several physiological processes. Recently, TMEM16A/ 
anoctaminl and TMEM16B/ anoctamin2 have been shown to function as CaCCs, but very little information is avail- 
able on the structure-function relations of these channels. TMEM16B is expressed in the cilia of olfactory sensory 
neurons, in microvilli of vomeronasal sensory neurons, and in the synaptic terminals of retinal photoreceptors. 
Here, we have performed the first site-directed mutagenesis study on TMEM16B to understand the molecular 
mechanisms of voltage and Ca^^ dependence. We have mutated amino acids in the first putative intracellular 
loop and measured the properties of the wild-type and mutant TMEM16B channels expressed in HEK 293T cells 
using the whole cell voltage-clamp technique in the presence of various intracellular Ca^* concentrations. We mu- 
tated E367 into glutamine or deleted the five consecutive glutamates sggEEEEEsgo and 39gEYE4oi. The EYE deletion 
did not significantly modify the apparent Ca"^ dependence nor the voltage dependence of channel activation. 
E367Qand deletion of the five glutamates did not greatly affect the apparent Ca^^ affinity but modified the voltage 
dependence, shifting the conductance-voltage relations toward more positive voltages. These findings indicate 
that glutamates E367 and sggEEEEEjgo in the first intracellular putative loop play an important role in the voltage 
dependence of TMEM16B, thus providing an initial structure-function study for this channel. 



INTRODUCTION 

Ca^^-activated CP channels (CaCCs) are expressed in 
many cell types, where they play various physiological 
roles. For example, CaCCs are involved in fast block of 
polyspermy in Xenopus laevis oocytes, in the regulation 
of smooth muscle contraction, in fluid secretion from 
exocrine glands, in the control of excitability in cardiac 
myocytes, as well as in olfactory, taste, and phototrans- 
duction (Frings et al., 2000; Hartzell et al., 2005; Leblanc 
et al., 2005; Petersen, 2005; Wray et al., 2005; Bers, 2008; 
Kleene, 2008; Lalonde et al., 2008; Petersen and Tepikin, 
2008; Duran et al., 2010; Kunzelmann et al., 2011a). 

Despite the fact that CaCCs are broadly present in sev- 
eral tissues, their molecular identity had remained elusive 
until 2008, when three independent studies reported 
that the expression of TMEM16A/ anoctaminl was 
associated with CaCCs (Caputo et al., 2008; Schroeder 
et al., 2008; Yang et al., 2008). The TMEM 16 family com- 
prises 10 members, and another member of the family, 
TMEM16B/anoctamin2, has also been shown to func- 
tion as a CaCC when heterologously expressed in axolotl 
oocytes (Schroeder et al., 2008) or in HEK 293T cells 
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(Pifferi et al., 2009; Stephan et al., 2009; Stohr et al., 
2009; Rasche et al., 2010; Sagheddu et al., 2010). 

The study of knockout mice for TMEM16A (Rock 
and Harfe, 2008) and for TMEM16B (Billig et al., 2011) 
further confirmed that CaCC activity was reduced or 
abolished in several cells (Flores et al., 2009; Galietta, 
2009; Hartzell et al., 2009; Kunzelmann et al., 2011b, 
2012; Huang et al., 2012; Pifferi et al., 2012; Sanders 
et al., 2012; Scudieri et al., 2012). 

Hydropathy analysis indicates that TMEM 16 proteins 
have eight putative transmembrane domains with both 
N- and C-terminal domains located at the intracellular side 
of the membrane, and the predicted topology has been 
experimentally confirmed for TMEM16G/anoctamin7 
(Dasetal., 2008). At present, TMEM 16A and TMEM 16B 
have been shown to function as CaCCs, whereas it is 
unclear whether the other members of the family are 
CaCCs (Duran and Hartzell, 2011; Huang et al., 2012; 
Scudieri et al., 2012). Furthermore, splice variants have 
been identified both for TMEM16A (Caputo et al., 
2008; Ferrera et al., 2009, 2011) and for TMEM16B 
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(Stepham et al., 2009). However, although the functional 
properties of different isoforms have been extensively 
investigated for TMEM16A, only preliminary data have 
been presented for TMEM16B (Saidu, S.P., A.B. Stephan, 
S.M. Caraballo, H. Zhao, and J. Reisert. 2010. Association 
for Chemoreception Sciences Meeting. Abstr. P68) . 

At present, very little is known about the structure- 
fiinction relations for these channels. The analysis of the 
sequence of TMEM 1 6A and TMEM16B did not reveal 
any canonical voltage-sensing or Ca^^-binding domains 
(Yang et al., 2008) , but a comparison among the biophys- 
ical properties of the TMEM16A splice variants pointed 
to the functional relevance of the first putative intracellu- 
lar loop (Caputo et al., 2008; Ferrera et al., 2009, 201 1) . 
Moreover, a recent study performed site-directed mu- 
tagenesis experiments on TMEM16A modifying some 
amino acids in the first putative intracellular loop and 
found that deletion of EAVK affected both the Ca^"^ and 
voltage dependence of TMEM16A (Xiao etal., 2011). 

Here, we aimed to perform a first site-directed muta- 
genesis investigation of TMEM16B to contribute to the 
understanding of the molecular mechanisms underly- 
ing the channel voltage and Ca^^ dependence. We iden- 
tified some acidic amino acids in the first intracellular 
loop of TMEM16B (sevE, sgeEEEEEsso, 399EYE401), which 
are conserved in TMEM16A, where some of them have 
been studied (Xiao et al., 2011). We mutated or deleted 
the indicated glutamates and made a comparison be- 
tween the electrophysiological properties measured 
in the whole cell configuration of the wild-type (WT) 
TMEM16B and its mutants. We have found that 367E 
and 386EEEEE390 contribute to the voltage-dependent 
regulation of the TMEM16B channel. 



MATERIALS AND METHODS 

Site-directed mutagenesis of TI\/IEIVI16B and 
heterologous expression 

Full-length mouse TMEM16B cDNA in pCMV-Sport6 mamma- 
lian expression plasmid was obtained from RZPD (clone identifi- 
cation, IRAVp968H1167D; NCBI Protein database accession no. 
NP_705817.1). Mutations were made using a PCR-based site- 
directed mutagenesis kit (Gene Tailor; Invitrogen) and con- 
firmed by DNA sequencing. HEK 29.ST cells (American Type 
Culture Collection) were transfected with 2 ]ig TMEM16B by 
using transfection reagent (FuGENE 6; Roche). Cells were co- 
transfected with 0.2 pg enhanced green fluorescent protein 
(eGFP; Takara Bio Inc.) for fluorescent identification of trans- 
fected cells. After 24 h, transfected cells were replated at a lower 
density and used for patch-clamp experiments between 48 and 
72 h from transfection. 

Electrophysiological recordings and ionic solutions 

Current recordings from HEK 293T cells expressing TMEM16B 
or its mutants were performed in the whole cell voltage-clamp 
configuration, as described previously (Pifferi et al., 2006, 2009). 
Patch pipettes were made of borosilicate glass (World Precision 
Instruments, Inc.) and pulled with a PP-830 puller (Narishige). 
Patch pipettes filled with the intracellular solution had a resistance 



of ^3-5 MO when immersed in the bath solution. Currents were 
recorded with an Axopatch ID or Axopatch 200B amplifier con- 
trolled by Clampex 9 or 10 via a Digidata 1332A or 1440 (Molecular 
Devices) . Data were low-pass filtered at 5 kHz and sampled at 10 kHz. 
Experiments were performed at room temperature (20-25°C). 
As reported previously (Pifferi etal., 2006), control experiments in 
nontransfected and only eGFP-transfected cells did not show any 
significant Ca^*-activated current. 

The standard extracellular solution contained (in mM): 140 
NaCl, .5 KCl, 2 CaClj, 1 MgClj, 10 glucose, and 10 HEPES, adjusted 
to pH 7.4 with NaOH. The intracellular solution filling the patch 
pipette contained (in mM): 140 CsCl, 10 HEPES, and 10 HEDTA, 
adjusted to pH 7.2 with CsOH, and no added Ca^* for the nomi- 
nally 0 Ca^* solution, or various added Ca^"^ concentrations, as 
calculated with the program WinMAXC (Patton et al., 2004), to 
obtain free Ca^* in the range between 0.5 and 100 pM. The free 
Ca^* concentrations were experimentally determined by Fura-4F 
(Invitrogen) measurements by using a luminescence spectropho- 
tometer (LS-50B; PerkinElmer), as described previously (Pifferi 
et al., 2006). The total CP concentration was 1.58 mM in the extra- 
cellular solution, whereas in the pipette solution it ranged from 
140 mM in 0 Ca^* to 160 mM in 100 pM Ca^+, with a calculated 
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Figure 1. I-V relations of TMEM16B. (A) Representative whole 
cell voltage-clamp recordings obtained with an intracellular solu- 
tion containing nominally 0 Ca^* or 1.5 pM Ca^*, as indicated. 
Voltage steps of 200-ms duration were given from a holding volt- 
age of 0 mV to voltages between — 100 and +100 mV in 20-mV 
steps, followed by a step to — 100 mV, as indicated in the top part 
of the panel. (B) Steady-state I-V relation measured at the end 
of the voltage steps (circles) or instantaneous I-V measured at 
the beginning of each voltage step (inverted triangles) from the 
cell shown in B. (C) Representative recordings at 1.5 pM Ca^* ob- 
tained with a voltage protocol consisting of a prepulse to +100 mV 
from a holding voltage of 0 mV, followed by voltage steps between 
— 100 and +100 mV in 20-mV steps, as shown in the top part of 
the panel. (D) I-V relations measured from tail currents (inverted 
triangles) or at the steady state (circles) . 
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equilibrium potential for CP of — 1.5 and +1.9 mV, respectively. 
All chemicals, unless otherwise stated, were purchased from 
Sigma-Aldrich. 

In most experiments, we applied voltage steps of 200-ms duration 
from a holding potential of 0 mV ranging from — 1 00 to + 1 00 mV 
(or from -200 to +200 mV), followed' by a step to -100 mV. 
A single-exponential function was fitted to tail currents to extrapo- 
late the current value at the beginning of the step to —100 mV. 
In another set of experiments, channels were activated by a 200-ms 
pulse to +100 mV, and then rapidly closed by the application of 
hyperpolarizing steps. Single-exponential functions were fitted to 
tail currents at each voltage step. 

Membrane capacitance and series resistance were compen- 
sated with the amplifier during the experiments. Membrane cur- 
rent density was calculated by dividing the current by the cell 
capacitance. The conductance, G, was calculated as G= // (V — V„v), 
where lis the tail current, V is the membrane voltage, and V,ev is 
the current reversal potential. Because in our experimental con- 
ditions the calculated equilibrium potential for CP ranged be- 
tween — 1.5 and +1.9 mV and the measured Vrev was close to 0 mV, 
Vrev was set to 0 mV in all calculations. 

Data analysis 

Data are presented as mean ± SEM, with n indicating the number 
of cells. Statistical significance was determined using paired or 
unpaired i tests or ANOVA, as appropriate. When a statistically 
significant difference was determined with ANOVA, a post-hoc 
Tukey test was done to evaluate which data groups showed 
significant differences. P-values of <0.05 were considered signi- 
ficant. Data analysis and figures were made with Igor Pro soft- 
ware (WaveMetrics). 
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Figure 2. Ca^*-dependent rectification of TMEMI6B. (A) Whole 
cell currents activated by the indicated [Ca^^Jj. Voltage protocol 
as in Fig. 1 A. (B) Average steady-state I-V relations from several 
cells {n = 3-6). (C) Average ratios between steady-state currents 
measured at +100 and —100 mV at various [Ca^*]i ( n = 3-6). 



RESULTS 

TMEM16B activation by Ca^* and voltage 
To study the activation of TMEM16B by [Ca^^]; and 
voltage, we performed whole cell voltage-clamp record- 
ings on HEK 293T cells transiently transfected with 
TMEM16B using intracellular solutions containing dif- 
ferent free [Ca^^Ji. Fig. 1 A shows that voltage steps be- 
tween — 100 and +100 mV from a holding voltage of 
0 mV elicited very small currents with a nominally 0-Ca^^ 
pipette solution (8 ± 3 pA/pF at +100 mV; n = 8) , whereas 
it induced large outward currents in the presence of 
1.5 pM Ca'\ 

In the presence of Ca^^, depolarizing voltage steps 
elicited an instantaneous outward current, indicating that 
channels were open at the holding potential of 0 mV, fol- 
lowed by a time-dependent outward relaixation (see also 
Fig. 5). Hyperpolarizing voltage steps induced instan- 
taneous inward currents followed by a relaxation toward 
less negative values, in agreement with previous results 
(Pifferi et al., 2009; Stohr et al., 2009; Rasche et al., 
2010). The I-V relation measured at the steady state 
showed a pronounced outward rectification, whereas 
the instantaneous I-V curve measured at the beginning 
of each step was linear (Fig. 1 B) . A similar result was 
obtained by activating TMEM16B with a different volt- 
age protocol: channels were first activated by a 200-ms 
prepulse to +100 mV, and then tail currents were in- 
duced by voltage steps between —100 and +100 mV in 
20-mV steps (Fig. 1 C). The I-V relation obtained by 
plotting the tail currents measured at the beginning of 
each step versus the step voltage was linear, whereas the 
steady-state I-V curve showed an outward rectification 
(Fig. 1 D) , as in Fig. 1 B. These results clearly demonstrate 
that the I-V relation of the open channel is linear, and 
therefore the outward rectification is a result of a voltage- 
dependent mechanism that favors channel opening at 
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Figure 3. Ca^* sensitivity of TMEMI6B. (A) Conductance den- 
sity calculated from tail currents measured at —100 mV after pre- 
pulses between — 100 and +100 mV as indicated was plotted versus 
[Ca^*]i (n = 3-6). Voltage protocol as in Fig. 1 A. Lines are the fit 
to the Hill equation (Eq. 1). (B) K1/2 and nn (inset) values plotted 
versus voltage. 
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depolarizing voltages. Thus, TMEM16B is activated by 
[Ca^*]i and modulated by voltage at low [Ca^^] i. 

To further examine the interplay between [Ca'^Ji and 
voltage in channel activation, we varied [Ca^^Ji (Fig. 2 A). 
Steady-state I-V relations meaisured at low [Ca^^Ji showed 
an outward rectification that became less pronounced 
as [Ca^^li increased (Fig. 2 B) . We calculated a rectifica- 
tion index as the ratio between the steady-state current 
at +100 and -100 mV at each [Ca^^^Ji. The rectification 
index was 4.8 ± 0.2 at 1.5 pM Ca^^ and decreased to 1.4 ± 
0.2 at 100 \iM Ca.-*, showing that the I-V relation is 
Ca^* dependent and becomes more linear as [Ca^^Ji 
increases (Fig. 2 C). 

To analyze the Ca^^ dependence of TMEM16B activa- 
tion at various voltages, we measured the dose-response 
relations. Tail currents at each [Ca'^Ji were measured at 
the beginning of the step to — 1 00 mV after prepulses 
ranging from ~ 100 to +100 mV. Fig. 3 A shows the aver- 
age conductance densities plotted versus [Ca'^Ji and fit 
at each voltage by the Hill equation: 

G = G_[Ca-];"/([Ca-X%^,/-), (D 

where G is the current density, Gmax is the maximal cur- 
rent density, K1/2 is the half-maximal [Ca"*]i, and uh is 
the Hill coefficient. 

The Hill coefficient was not voltage dependent with a 
value of 1.2 at -100 mV and 1.1 at +100 mV. The finding 
that the Hill coefficient was >1 indicates that the bind- 
ing of more than one Ca^* ion is necessary to open the 
channel. Ki/<i slighdy decreased with membrane depo- 
larization from 4.3 pM at - 1 00 mV to 1 .6 pM at +100 mV, 
as illustrated in Fig. 3 B. These data show that the Ca^* 
sensitivity of TMEM16B is moderately voltage dependent, 
in agreement with previous results obtained with inside- 
out patches (Pifferi et al., 2009; Stephan et al., 2009). 
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Figure 4. Voltage dependence of TMEM16B. (A) Normalized 
conductances at the indicated [Ca^^]; calculated from tail cur- 
rents at —100 mV after prepulses between —200 and +200 mV 
were plotted versus the prepulse voltage (n = 4-9). Lines are the 
fit to the Boltzmann equation (Eq. 2) . (B) V1/2 and z (inset) values 
plotted versus [Ca^*]i. 



The voltage dependence of steady-state activation 
(G-V relation) was analyzed by measuring tail currents at 
the beginning of a step to —100 mV after prepulse volt- 
ages between —200 and +200 mV. The range of voltages 
was extended from the previous voltage protocols to ob- 
tain a better estimate of voltage dependence. Fig. 4 A 
shows the average conductance activated at a given 
[Ca^^li plotted versus membrane voltage and fit by the 
Boltzmann equation: 

G/G_= l/{l + exp[z{V,,,-V)F/RT]}, (2) 

where G/ G^ax is the normalized conductance, z is the equiv- 
alent gating charge associated with voltage-dependent 
channel opening. Vis the membrane potential, V1/2 is 
the membrane potential producing half-maximal acti- 
vation, F is the Faraday constant, R is the gas constant, 
and T is the absolute temperature. 

The maximal conductance density G^ax was determined 
by a global fit of G-V relations, and G at each [ Ca^*] i was 
then normalized to the same Gmax- Because at the smaller 
[Ca^^li the prediction of G^ax from the fit could be 
affected by a large error, we also estimated G^ax at each 
[Ca^+li. G„ax at 0.5 pM Ca^^ was 4.1 ± 0.4 nS/pF, not 
significandy different from the value of 4.7 ± 0.4 nS/ pF 
at 100 pM Ca^*, indicating that the estimate of G^ax was 
litde affected by [Ca^*];. Fig. 4 A shows that increasing 
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Figure 5. Activation and deactivation kinetics of TMEM16B. 
(A) Representative recordings at the indicated [Ca^*];. Voltage 
protocol as in Fig. 1 A, with voltage steps from a holding voltage of 
0 between +40 to +100 mV in 20-mV steps. Red dashed lines are 
the fit to a single-exponential function. (B) Average activation time 
constants (T^ct) plotted versus voltage (ra = 6-8). (C) Represen- 
tative recordings at the indicated [Ca^*]i. Voltage protocol as in 
Fig. 1 C, with a prepulse to +100 mV and tail currents indviced by 
voltage steps between —100 and +100 mV in 20-mV steps. Only 
tail currents are illustrated. Red dashed lines are the fit to a single- 
exponential function. (D) Average deactivation time constants 
(Tdeact) plotted versus voltage (m=4— 9). 



288 Voltage dependence of TMEM16B/anoctamin2 



[Ca^^li produced a leftward shift in the G-V relation: 
Vi/2 was 124 ± 20 mV at 1.5 \iM Ca?'' and became -115 ± 
18 mV at 100 pM Ca^*, whereas the equivalent gating 
charge was not largely modified (z = 0.23-0.30). Thus, 
Vi/2 decreased as [Ca^^Ji increased, indicating that 
more channels can be activated by depolarization in 
the presence of a high [Ca^^Ji (Fig. 4 B). At a given 
[Ca^^li, the conductance increased with depolarization, 
showing that the conductance depends both on [Ca^^Ji 
and voltage. 



Activation and deactivation kinetics are regulated by 
[Ca^*]i and voltage 

To characterize activation and deactivation kinetics, we 
analyzed the time-dependent components in response to 
voltage steps in the presence of a given [Ca^^Ji. As shown 
in Figs. 2 A and 5 A, current activation in response to 
depolarizing voltage steps had two components: an in- 
stantaneous time-independent current, related to the 
fraction of channels open at the holding voltage of 0 mV, 
followed by an outward time-dependent relaxation, a 
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Figure 6. TMEM16B mutations. (A) Predicted topology of TMEM16A and TMEM16B from hydropathy analysis. (B) Alignment be- 
tween mouse TMEM16A {a,c, available from GenBank/EMBL/DDBJ under accession no. NM_178642.4) and the retinal isoform of 
TMEM16B used in this study (NP_70.5817.1), with the mutations or deletions highlighted in color. (C-F) Representative recordings at 
the indicated [Ca^*]i for E.S67Q (C), AE5 (D), AEYE (E) mutants, and WT (F). Traces for WT are the same as in Fig. 2 A. Voltage proto- 
col as in Fig. 1 A. (G) I-V steady-state relations (n = 3-8). (H) Average ratios between currents measured at +100 and —100 mV plotted 
versus [Ca^*]i for each mutant (n = S-8). 
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result of the increaise in the fraction of channels opened 
by depolarization. The time-independent component 
became larger as voltage or [Ca"^]i increaised. 

To examine the activation kinetics, we analyzed the 
time-dependent component of the current elicited by 
depolarizing voltage steps. Fig. 5 A shows that most of 
the time course of time-dependent relaxations was 
well fit by a single-exponential function. The time 
constant of current activation, Tact, in the presence of 
1.5 pM Ca'+ was 8.1 ± 0.8 ms at +100 mV and did not 
vary as a function of voltage at a given [Ca^*] i (Fig. 5 B) . 
At +100 mV, Tact at 100 pM Ca^" was 3.9 ± 1.4 ms, sig- 
nificantly smaller than the value of 8.1 ± 0.8 ms at 
1.5 pM Ca^*, showing that an increase in [Ca^^Ji accel- 
erated activation. 

The time constant of current deactivation (Tdeact) 
was calculated by fitting with a single-exponential 
function the tail currents obtained after a prepulse at 
+100 mV by voltage steps ranging between —100 and 
-20 mV (Fig. 5 C). In the presence of 0.5 pM Ca^^ 
Tdeact was 3.0 ± 0.2 ms at -100 mV and 5.4 ± 0.5 ms at 
— 20 mV, showing that less negative voltages slowed 
deactivation (Fig. 5 D). At -100 mV, Tjeact at 1.5 pM 
Ca^^ was 7.2 ± 0.8 ms, significantly different from the 
value of 3.0 ± 0.2 ms at 0.5 pM Ca^^, showing that an 
increase in [Ca^*]i slowed deactivation. 
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Figure 7. Ca'* sensitivity of TMEM16B mutants. Conductance 
density calculated from tail currents measured at — 1 00 mV after 
prepulses between —100 and +100 mV as indicated was plotted 
versus [Ca^*]; for E367Q (A; n = 3-6), AE5 (B; n = 3-5), and AEYE 
(C; n = 3-8) mutants. Lines are the fit to the Hill equation (Eq. 1 ) . 
(D) Ki/2 values plotted versus voltage for each mutant. 



In summary, the activation kinetics are voltage in- 
dependent and become faster by increasing [Ca^^Ji, 
whereas the deactivation kinetics are prolonged by de- 
polarization and by increasing [Ca'*]i. 

Functional characterization of mutations in the first 
putative intracellular loop 

To investigate the molecular mechanisms responsible for 
channel activation by Ca'* and by voltage, we performed 
a site-directed mutagenesis study. Hydropathy analysis in- 
dicates that each member of the TMEM16 family has 
eight transmembrane domains (Fig. 6 A) . Analysis of the 
sequence of TMEM16B does not reveal the presence of 
any typical voltage sensor or Ca^^-binding domain. How- 
ever, some acidic amino acids are located in the first puta- 
tive intracellular loop between transmembrane segment 
2 and 3, and we hypothesized that some of them may be 
involved in Ca^* and/ or voltage activation of TMEM16B. 
As illustrated in Fig. 6 B, we mutated glutamate at position 
367 into glutamine (E367Q) , deleted the five consecutive 
glutamate residues sgeEEEEEago ( AE5) , or deleted 399EYE401 
(AEYE), and meaisured their biophysical properties. 

Fig. 6 (C-F) illustrates recordings from each mutant 
channel in the presence of various [Ca'^Jj. Similar to WT 
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Figure 8. Voltage dependence of TMEM16B mutants. Normal- 
ized conductances at the indicated [Ca^*]i calculated from tail 
currents at — 100 mV after prepulses between —200 and +200 mV 
were plotted versus the prepulse voltage. Black lines are the fit 
to the Boltzmann equation (Eq. 2) for WT from Fig. 4 at 100 pM 
Ca^* (solid line) or at 1.5 pM Ca^* (dashed line) . Colored lines are 
the fits to the Boltzmann equation for E367Q (A; n = 3-4), AE5 
(B; n = 3-5), and AEYE (C; n = 3-6) mutants. (D) Average 
values plotted versus [Ca^*]i. 
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(Fig. 2 A), the steady-state I-V relation for each mutant 
was Ca'* dependent, showing an outward rectification 
at low [Ca^^li that became less pronounced as [Ca^*]i 
increased (Fig. 6 G). However, although the overall 
Ca^* dependence was similar, the rectification index, 
measured from the ratio between steady-state currents 
at +100 and —100 mV, was significantly higher at every 
[Ca^^li in E367Qand AEg mutants than in WT, whereas 
it remained similar in AEYE mutant channel (Fig. 6 H) . 

The dose-response relations for each mutant chan- 
nel, evaluated from tail currents as described previously 
for the WT channel (Fig. 3) , were fit by the Hill equa- 
tion (Fig. 7, A-C). Fig. 7 D shows that at +100 mV 
( - 1 00 mV) was 1 .6 \M (4.3 pM) in WT, 1 .3 pM (4.0 pM) in 
E367Q, 2.2 pM (4.0 pM) in AE5, and 1.3 pM (3.2 pM) in 
AEYE. The Hill coefficient uh at +100 mV (~100 mV) 
was 1.1 (1.2) inWT, 1.6 (1.2) inE367Q, 1.4 (2.9) in AE5, 
and 2.0 (1.7) in AEYE. Thus, the mutations produced only 
some very small changes in Ki/^ or Uh, but overall no 
strong modifications in Ca"* sensitivity were observed. 

The G-V relations in mutant channels were measured 
at each [Ca'^Ji and compared with the corresponding 
relations in WT channels. Fig. 8 A shows that the E367Q 
mutation produced a rightward shift of the G-V relation 
at a given [Ca^^Ji with respect to WT; indeed, V1/2 
changed from 124 + 20 mV in WT to 169 ± 6 mV in 
E367Q at 1.5 pM Ca^", and from ~115 ± 18 mV in WT 
to 44 ± 8 mV in E367Q at 100 pM Ca'+ (Fig. 8 D) . The 
deletion AE5 also shifted the G-V relations to the right 
(Fig. 8, B and D) : Vi/, changed from 124 + 20 mV in WT 
to 248 + 39 mV in AE5 at 1.5 pM Ca.^\ and from -115 + 
18 mV in WT to 58 + 15 mV in AE5 at 100 pM Ca^^ Dif- 
ferendy from the previous mutants, the AEYE deletion 
did not produce any significant change in the G-V rela- 
tions (Fig. 8, C and D) . The equivalent gating charge for 
each mutant varied between 0.15 and 0.32, values similar 
to those of the WT channel (z= 0.23-0.30). Thus, E367Q 
and the AE5 deletion modified the voltage sensitivity: at a 
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Figure 9. Activation and deactivation kinetics of TMEM16B mu- 
tants. Kinetics were measured as explained in Fig. 5. (A) Average 
activation time constants (Tact) plotted versus voltage for E367Q 
(n = 5), AE5 (b= 3), and AEYE (n = 6) mutants. (B) Average deac- 
tivation time constants (T^eact) plotted versus voltage for E367Q (n = 
4), AE5 (ra = 4), and AEYE (k = 5) mutants. 



given [Ca^^li, fewer channels can be open by depolariza- 
tion compared with WT. 

The kinetic properties of activation and deactivation of 
mutant channels also showed some interesting changes 
compared with WT channels. Upon depolarizing volt- 
age steps, the activation of mutant channels was still 
characterized by two components: an instantaneous time- 
independent current, followed by an outward time- 
dependent relaxation (Fig. 6), which was well fit by a 
single-exponential function as in WT channels. In the 
presence of 1.5 pM Ca^^, T^ct at +100 mV was 2.8 + 0.3 ms 
in E367Q, faster than 7.5 + 0.7 ms in the WT channel, 
whereas it became slower than WT in AE5 (17.7 + 3.0 ms) 
and in AEYE (25.5 + 2.3 ms). These results indicate that 
each mutation altered the time course of activation. 
Indeed, the time necessary to respond to a depolarization 
decreased in E367Q, whereas it was progressively pro- 
longed in AE5 and in AEYE compared with WT. As in the 
WT channel, r^a in each mutant was not significantly 
modified by voltage (Fig. 9 A) . 

Deactivation kinetics was also well fit by a single- 
exponential function and, similarly to WT, Tjeact showed 
an increase at less negative voltages for each mutant 
channel (Fig. 9 B). In the presence of 1.5 pM Ca^*, Tdeact 
at -100 mV was 1.6 + 0.3 ms in E367Q, smaller than 7.2 + 
0.8 ms in the WT channel, whereais it was not significantly 
differentfrom WT in AE5 (6.8 + 0.3 ms) and became larger 
than WT in AEYE (12.3 +1.5 ms) . The time necessary for 
channels to close upon repolarization decreaised in E367Q 
but remained similar in AE5, and it wais prolonged in AEYE 
compared with WT. Thus, E367Q and AEYE mutants also 
showed a modified time course of deactivation. 



DISCUSSION 

Here, we have provided the first site-directed mutagen- 
esis study to investigate structure-function relations of 
the TMEM16B channel. Because previous studies have 
shown that TMEM16B in excised inside-out patches has 
a significant nmdown (Piffeti et al., 2009, Fig. 5; Stephan 
et al., 2009, Fig. 3 A) , whereas whole cell recordings are 
rather stable (Pifferi etal., 2009, Fig. 1 h), we decided to 
use the whole cell configuration. 

We first characterized the WT TMEM16B channel and 
established one important difference between TMEM16A 
and TMEM16B activation properties in the absence of 
[Ca^^^li. Indeed, we found that TMEM16B cannot be acti- 
vated by voltages up to +200 mV in the absence of Ca^* 
(32 +10 pA/ pF; w = 6; not depicted) , whereais recent data 
from HartzeU's laboratory showed that TMEM16A was ac- 
tivated by strong depolarization in the absence of Ca^^ 
(-^140 pA/pF at +200 mV; Fig. 5 A in Xiao et al., 2011). 
Thus, our data show that TMEM16B needs Ca'^ to be acti- 
vated differently from TMEM16A, which can be activated 
by voltage also in the absence of Ca^^ (Xiao etal., 2011). 
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In the presence of Ca~*, dose-response relations for 
TMEM16A and TMEM16B obtained by different labo- 
ratories reported variable values for Ki/^. For TMEM 1 6A, 
from inside-out recordings, K1/2 at +60 mV ( — 60 mV) 
wasO-SpM (2.6 pM) (Yang etal., 2008), and at + 100 mV 
(-100 mV) it was 0.4 pM (5.9 pM) (Xiao et al., 2011), 
whereais from whole cell recordings at+100 mV ( — 40 mV) 
it was 332 nM ('-TOO nM) (Ferrera et al., 2009). For 
TMEM16B, from previous work in inside-out patches, K1/2 
at +50 mV (-50 mV) was 3.3 pM (4.9 pM) (Pifferi et al., 
2009), and at+40mV (-40 mV) it was 1.2 pM (1.8 pM) 
(Stephan et al., 2009), whereas from whole cell record- 
ings we found that K1/2 at +40 mV (-40 mV) was 2.0 pM 
(2.7 pM), and at+100 mV (-100 mV) it was 1.6 pM 
(4.3 pM) (Fig. 3). Although there are some differences 
among studies reported from different laboratories, 
every report showed that the apparent affinity for Ca"* 
is slightly voltage dependent, with higher apparent Ca"^ 
affinity at positive voltages, and the Hill coefficients are 
consistently higher than one, indicating that more than 
a Ca^^ ion is necessary to activate the channels. A compari- 
son between TMEM16A and TMEM16B shows a fourfold 
difference between K1/2 values at +100 mV: 0.4 pM (Xiao 
et al., 2011) for TMEM16A and 1.6 pM for TMEM16B 
(Fig. 3), indicating a lower apparent affinity for Ca'* of 
TMEM16B compared with TMEM16A. 

A critical question about the function of TMEM16A 
and TMEM16B is the following: what are the molecular 
mechanisms responsible for Ca^^ and voltage modula- 
tion of channel gating in each channel? Galietta's labo- 
ratory (Ferrera et al., 2009) has shown that human 
TMEM16A has various protein isoforms generated by 
alternative splicing, and it has labeled the four identi- 
fied alternative segments as a, b, c, and d. A rare minimal 
version of TMEM16A lacking all alternative segments, 
TMEM16A (0), still shows CaCC properties, although the 
voltage dependence is reduced, (Caputo et al., 2008; 
Ferrera et al., 2009, 2011). Ferrera et al. (2009) showed 
that segment b modified the Ca'^ sensitivity by nearly four- 
fold, decreaising the apparent half-effective concentration 
at +80 mV fi-om 350 to 90 nM, whereas segment c affected 
the voltage dependence but not the Ca^^ sensitivity of 
human TMEM16A {abc). Segment c is composed of the 
four amino acids EAVK, which have also been recendy 
deleted from mouse TMEM16A {ac) in a study from 
Hartzell's laboratory (Xiao et al., 2011). Differently from 
Ferrera et al. (2009), Xiao et al. (2011) found that dele- 
tion of EAVK modified both Ca^^ and voltage dependence 
of TMEM16A. The discrepancy between the results can be 
a result of differences between human TMEM16A {ab) 
and mouse TMEM16A (a), and/ or to the different tech- 
niques, whole cell versus inside-out recordings, used for 
the experiments in the different laboratories. Although 
the two studies reached some different conclusions, they 
both pointed to the relevance of the segment cin the regu- 
lation of the TMEM 16A functional activity. 



TMEM16B is expressed in the retina, at the synaptic 
terminal of photoreceptors (Stohr et al., 2009; BUlig et al., 
2011), in the cilia of olfactory sensory neurons, and in 
the microvilli of vomeronasal sensory neurons (Stephan 
et al., 2009; Rasche et al., 2010; Sagheddu et al., 2010; 
Billig et al., 2011; Pifferi et al., 2012). Zhao's laboratory 
showed that the major TMEM16B olfactory isoform dif- 
fers from the retinal isoform in the absence of the exon 
encoding the four amino acids ERSQ in the first puta- 
tive intracellular loop (Stephan et al., 2009). It is worth 
pointing out here that segment c (EAVK.) in TMEM16A 
is not present in TMEM16B, but that ERSQ residues are 
located in the corresponding positions in the retinal 
isoform of TMEM16B (Fig. 6). A comparison between the 
biophysical properties measured in inside-out patches 
from the retinal isoform (Pifferi et al., 2009) and from 
the olfactory isoform (missing ERSQ Stephan etal., 2009) 
did not reveal any major difference in the rectification 
properties and in the dose-response relations between 
the two isoforms, although we cannot exclude that 
more detailed biophysical studies may reveal subtle dif- 
ferences. Indeed, the functional properties of addi- 
tional isoforms for TMEM16B are under investigation 
(Saidu, S.P., AB. Stephan, S.M. Caraballo, H. Zhao, and 
J. Reisert. 2010. Association for Chemoreception Sciences 
Meeting. Abstr. P68) . 

Although the amino acidic sequences of both TMEM16A 
and TMEM16B lack any classical voltage-sensor or Ca^*- 
binding domain, a series of five consecutive glutamates 
located in the first putative intracellular loop has been 
identified as a good candidate to play a role in channel 
gating. Moreover, we have investigated if other glutamates 
in the same loop could also be involved in the activation 
of TMEM16B by Ca^"^ and voltage. We found that dele- 
tion of the five glutamates, AE5, did not gready affect 
the apparent affinity for Cd?* (Fig. 7) , but it significantly 
shifted the activation curve to the right. Indeed, V1/2 at 
1.5 pM Ca^* changed from 124 mV in WT to 248 mV, 
whereas the equivalent gating charge was not modified. 
In addition, the time necessary to respond to a depolar- 
ization was prolonged in AE5, whereas the deactivation 
constant was not significantly affected (Fig. 9) . Thus, the 
five consecutive glutamates are involved in the voltage 
dependence of the TMEM16B channel, whereas they do 
not seem to play a significant role in the apparent affinity 
for Ca^^. These results are in agreement with a recent 
study in TMEM16A, showing that the substitution of the 
four correspondent glutamates into alanines (444EEEE/ 
AAAA447) did not gready affect the apparent affinity for 
Ca^* but modified the voltage dependence, producing a 
shift of the activation curve to the right (Xiao et al., 201 1 ) . 

In the TMEM16B mutant E367Q, both activation and 
deactivation kinetics were shortened; the dose-response 
relation for Ca^* was not strongly modified, while the 
activation curve was shifted to the right. Finally, the de- 
letion AEYE produced an increase in the time constants 
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for activation and deactivation, whereas it did not cause 
any large change in apparent affinity for Ca^^ or in volt- 
age sensitivity. 

Collectively, our results indicate that glutamates E367 
and 386EEEEE390 in the first putative intracellular loop 
play a relevant role in the modulation of the voltage 
dependence of TMEM16B. 

Conclusions 

In conclusion, we have found evidence that the five con- 
secutive glutamates in the first putative intracellular loop 
are not involved in Ca^^ sensitivity in TMEM16B but 
have an important role in voltage dependence. Another 
glutamate in position 367 plays a similar role, further 
indicating that the first intracellular loop is involved in 
voltage-dependent activation of TMEM16B. 

At present, the location of the Ca^^-binding site in 
TMEM16A and TMEM16B remains unknown. It is pos- 
sible that several residues in different regions contrib- 
ute to bind Ca^^ ions, but it cannot be excluded that the 
Ca^^-binding site is located in an accessory subunit ex- 
pressed both in HEK 293T cells and in axolotl oocytes. 
Future work will have to shed light on the intricate 
mechanisms that couple Ca^* gating and voltage depen- 
dence, including intriguing interactions between gating 
and permeation. 
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